Internal Extinction in Spiral Galaxies in the Near Infrared 
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ABSTRACT 

In order to study the effects of internal extinction in spiral galaxies we searcli for correlations 
of near infrared (NIR) photometric parameters with inclination. We use data from the 2 Micron 
All-Sky Survey (2MASS) Extended Source Catalog (XSC) on 15,224 spiral galaxies for which we 
also have redshifts. For 3035 of the galaxies, I-band photometry is available which is included in 
the analysis. From the simple dependence of reddening on inclination we derive a lower limit to 
the difference in magnitude between the face-on and edge-on aspect of 0.9, 0.3 and 0.1 magnitudes 
in I (0.81/im), J (1.25 /zm) and H (1.65 /im) bands. We find that the faintest isophotal radius 
reported in the XSC (at the 21st magarcsec"^ level) is closer to the centers of the galaxies 
than other common isophotal measures {e.g. the 23.5 magarcsec"^ radius in I-band), and argue 
that it should not be assumed to represent an outer isophote at which galaxies are transparent 
at all viewing angles. A simple linear extinction law (i.e. AM — j\og(a/b)) is not adequate 
for the full range of disk inclinations and we adopt both a bi-linear and a quadratic law. A 
simple photometric model is used to explain the observed behavior. Internal extinction depends 
on galaxy luminosity. We show that for galaxies with a Kg total magnitude dimmer than -20, 
-20.7 and -20.9 the data indicates zero extinction in J, H and Kg respectively, while disk opacity 
increases monotonically with increasing disk luminosity. 

Subject headings: galaxies: fundamental parameters - galaxies: photometry - galaxies: spiral - galaxies: 
statistics - dust, extinction - infrared radiation 



1. Introduction 

The determination of an accurate value for the 
internal extinction in galaxies is of relevance both 
in the understanding of their structure and in ap- 
plications such as the determination of redshift 
independent distances. Edge-on spirals are used 
preferentially in studies that apply the luminosity- 
linewidth relation (TuUy & Fisher 1977) and these 
are, of course, the galaxies in which extinction ef- 
fects will be the largest. Inadequate corrections 
for extinction will produce systematic distance er- 
rors in our catalogs. Recent discussions of this 
issue can be found in Disney & Burstein (1995) 
and Calzetti (2001), while a comprehensive his- 
tory of the study of internal extinction in galax- 
ies can be found in the introduction of Huizinga 
(1994). Analyses similar to the ones carried out in 
this paper can be found in Giovanelli et al. (1994; 



hereafter G94, and 1995) and TuUy et al. (1998). 

Since there is intrinsic variability in the proper- 
ties of galaxies, the determination of the amount 
of extinction is most readily obtained in a statisti- 
cal way, in which case large samples are desirable. 
However, sample selection effects need to be well 
understood in order to avoid interpreting spurious 
effects as astrophysically significant. There will be 
further discussion of the effects of such biases on 
the results of this paper in Sections 5.3 and 6.2. 

In this paper we report on statistical tests which 
search for the effects of extinction in spiral galax- 
ies in the near infra-red (NIR) bands, I (0.81 /xm), 
J (1.25 Aim) , H (1.65 ^nJl) and Kg (2.17 ^m). NIR 
bands are often used in cases where accurate total 
magnitudes of galaxies are needed, as they clearly 
minimize, if not eliminate, the need for extinction 
corrections. For this reason it is particularly in- 
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tercsting to search for the effects of extinction in 
these bands. 

Samples and available data are described in 
Sections 2 and 3. In Section 4 we apply simple 
tests for reddening with inclination and in Sections 
5 and 6 we discuss in more detail the effects of 
extinction on the isophotal radius and total mag- 
nitudes respectively. A simple photometric model 
for the extinction is discussed in Section 7 and in 
Section 8 luminosity dependences are discussed. A 
summary of the conclusions is given in Section 9. 

2. Description of the Samples 

The Two Micron All-Sky Survey (2MASS) com- 
pleted its coverage of the near infra-red sky in 
February 2001. Data were taken in three NIR 
bands, J, H and Kg using two identical 1.3m tele- 
scopes, one in Arizona, the other in Chile. With an 
extended source sensitivity of 14.8, 14.0 and 13.4 
mag at J, H and Kg (Jarrett et al. 2000) and an an- 
gular resolution of about 2-3 arcseconds, 2MASS 
has detected more than 1.65 million galaxies (Jar- 
rett et al. 2003); the resultant data are available 
in the Extended Source Catalog (XSC)"'^. 

Two samples were extracted from the 2MASS 
All-Sky release XSC. For statistical studies the 
samples should be as large as possible, but it is 
also important for them to have well defined se- 
lection criteria. The first sample (and by far the 
largest) consists of all galaxies with well defined 
spiral morphology and known redshift taken from 
the Arecibo General Catalog (a private compila- 
tion by the authors known as the AGC). Nearby, 
peculiar motions contribute significantly to ob- 
served recessional velocities so we further limit the 
AGC to cz > 3000 km s^^ in order to use the red- 
shift as a reliable distance indicator. This large 
sample is referred to as the AGC2;. The AGC is an 
all-sky compilation, but has varying completeness 
functions in different areas of the sky. A second, 
smaller sample consists of galaxies in the AGG^; 
for which I-band photometry is also available (as 
will be described below). This sample we call the 
SFI2 as it contains all of the SFI data (Giovanelli 
et al. 1997; Hayncs et al. 1999a, b) plus simi- 
lar but newer and as yet unpublished additions to 

Available through the online catalog 
query service of the IR Science Archive: 
http://irsa.ipac.caltech.edu/applications/Gator/ 



that. This sample is a homogeneous all-sky sample 
of mostly Sc galaxies. In summary: 

AGC2 Spiral galaxies of known redshift from the 
Arecibo General Catalog. The redshift was 
further confined to cz > 3000 km s~^. This 
sample includes a total of 15244 galaxies. 

SFI2 Galaxies from AGCz which also have I- 
band photometry available. This sample has 
a total of 3035 galaxies, which are mostly of 
types Sc and Sbc. 

The detection fractions of these galaxies in the 
2MASS XSC (within 12" of the previously cata- 
loged position) are 80% for the AGCz and 85% 
for the SFI2 sample. We expect that this is due 
in part to the detection limits of the 2MASS sur- 
vey and in part to the accuracy of the currently 
tabulated positions of the galaxies in the AGC. In 
support of this assumption the detection fraction 
increases for the SFI2 galaxies that are expected 
to have more accurate positions in the AGC. For 
the Second Incremental release of 2MASS data, a 
discussion of completeness is included in the on- 
line Explanatory Supplement^. When the search 
for a match is extended to a radius of 0.5' of 
known galaxies from the Zwicky Catalog, the de- 
tection rate is 95%. In the NIR, surface bright- 
ness decreases towards later type galaxies (Jarrett 
et al. 2003) making late type, gas rich galaxies 
(like the Sc galaxies which make up the bulk of 
the SFI2 sample) harder to detect than E/SO/Sa 
galaxies. 

3. Description of data 

For all of the galaxies in our samples vari- 
ous data are available. These included redshifts, 
coarse blue magnitudes and angular sizes from the 
AGC, photometry from 2MASS and I-band pho- 
tometry for the SFI2 sample. In the computa- 
tion of absolute magnitudes a value of Ho = lOO/i 
km s~^ Mpc~^ is adopted. 

3.1. I-band photometry 

The I-band photometry available for the SFI2 
sample is described in Haynes et al. (1999a). 

^Tlic Explanatory Supplement to the 2MASS Second Incre- 
mental Release (Curtri et al. 2000) is an on-line document: 
http://www.ipac.caltech.edu/2mass/releases/docs.html 
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Quantities of particular interest are: 

{a/b)i - I-band axial ratio derived from a fit of 

elliptical surface brightness contours to the 
I-band images. The average error in log(a/6) 
is 0.03. 

?"23.5 - radius measured at the I-band isophote of 
23.5 magarcsec"^. 

mi - total I-band magnitude. Calculated from an 
extrapolation of a fit to the disk to 8 scale 
lengths of the disk. Typical total errors are 
in the range 0.03-0.06 mag. 

rd - scale length of the disk derived from a fit 
to the one-dimensional I-band photometric 
profile. 

3.2. Photometric quantities from the 2MASS 

xsc 

The photometric quantities in the XSC are 
derived by inspecting the 2MASS images with 
the extended source processor (GALWORKS), 
which was developed for this purpose (see Jar- 
rett et al. 2000). By the time the 2MASS im- 
ages get to GALWORKS a point source spread 
function has already been derived for each image. 
GALWORKS identifies extended sources (galax- 
ies as well as other more nearby resolved sources) 
and measures a variety of parameters for each 
one as described in Jarrett et al. (2000) and in 
the 2MASS Explanatory Supplement (Curtri et 
al. 2000). Very large galaxies (greater than 5' in 
size at ~ 20th magarcsec"^ in J-band) are first 
identified in in a separate large galaxy atlas^ and 
then added to the XSC. The quantities from the 
XSC that we use are: 

(a/ 6) J - J-band axial ratio fit to the 3cr isophote 
(which is close to the 20th magarcsec"^ 
isophote). The center of this ellipse fit is 
fixed to the peak intensity pixel and since 
only one isophote is fitted, it is an approx- 
imation to the ellipticity of the galaxy. If 
the 3cr isophote is too close to the center of 
the galaxy the bulge might interfere with the 
derived ellipticity, making the galaxy appear 
more face-on than it really is. Of the three 

^See: http://irsa.ipac.caltech.edu/applications/2MASS/LGA/ 



2MASS bands, the J axial ratio should suf- 
fer the least from this effect; we note that 
no systematic offset is seen between the J- 
band and I-band axial ratios, as can be seen 
Figure 4d. An estimate of the average error 
of 0.1 in log(a/6) is derived from the scatter 
in this plot; this is more than three times as 
large as the average error estimated for the 
equivalent I-band quantity. 

r'li , r^i , r^i - semi-major axis of the isophotal el- 
liptical aperture set to the 21st magarcsec"^ 
isophote individually in each band. An 
alternative to this would have been the 
semi-major axis of the ellipse; fit to the 
20th magarcsec"^ isophote, which is the 
preferred aperture in the XSC (Jarrett et 
al. 2003). In Section 5.1 it is shown that 
r2i in J-band is at ~ 2.5 scale lengths 
of the disk (r^), while in H and Kg it is 
^ 3.3rd. A similar analysis for r2o shows 
that ~ 1.6r(i, and v^-^ ~ 2.2rd. We 
choose to use r2i in order to trace the disks 
as far out as possible (which still is not far 
compared to rjs.s ~ 3.5rd). Since the typ- 
ical 1 (J background noise is 21.4, 20.6 and 
20.0 mag in J, H and Kg respectively (Jar- 
rett et al. 2003), we must be careful in using 
the radii set to the 21st magarcsec"^ level, 
especially in H and Kg. In J-band r2i should 
be reasonably well defined, and does indeed 
show good correlation with r2o with a scatter 
of about 3 ". The same analysis for H and 
Kg shows a similar correlation with a slightly 
larger scatter (3.5" and 4 " respectively). 
When r is smaller than 7" or greater than 
120" (for either isophote and in all bands) it 
does not follow the correlation well. These 
radii are not used. 

J2if, H2if, K2if - magnitudes within the 21st 
magarcsec"^ elliptical isophote set in J- 
band (referred to as "J fiducial" are used 
to construct colors because of the necessity 
of having a consistent aperture. The un- 
certainty on these magnitudes listed in the 
XSC is that associated with Poisson noise 
and background subtraction and does not 
include other errors such as the zero point 
calibration and the uncertainty in the ellip- 
tical isophotal aperture (see Jarrett et al. 
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2000 for more details). This formal uncer- 
tainty is small (with the largest quoted at 
only 0.2 mag) and the real uncertainty is 
likely to be larger. The "Level One Sci- 
ence Requirements" of 2MASS (which can 
be found in the on-line Explanatory Supple- 
ment) require that the extended source pho- 
tometry be accurate and uniform to within 
10% , and we assume this is a closer estimate 
of the error in the magnitudes. 

J2if5 H2if, K2if - magnitudes measured within 
the 21st magarcsec"^ level set in J-band are 
used to construct colors. 

Jextt Hcxtj Koxt - XSC "extrapolated" magni- 
tudes designed to represent the total magni- 
tude of the galaxy. These magnitudes are 
derived from an extrapolation of the ra- 
dial surface brightness profile from the 20th 
magarcsec"^ elliptical isophote down to rtot 
which is set to 5 times the scale length from 
the J-band fit. See Jarrett et al. (2003) for 
more details. 

3.3. Corrections 

In order to use the photometry from the 2MASS 
XSC to its best potential for our purposes, various 
corrections for seeing, galactic extinction and cos- 
mological redshift need to be applied. These cor- 
rections have already been applied to the I-band 
quantities in a similar manner as that described 
below (for details, see Haynes et al. 1999a). 

3.3.1. Seeing correction 

Figure la shows a plot of log(a/6)j vs. the 
20th magarcsec^^ isophotal radius (the approx- 
imate isophotal radius at which a/6 is reported). 
We see that there is a lack of galaxies with both 
high axial ratio and small size. This is due mainly 
to the resolution of the images set by the seeing 
and the pixel size. 

In the XSC, seeing is characterized using a mod- 
ified exponential of the form 



/(r) = /(0)exp 



1//3' 



(1) 



"radial shape" (sh = q x /?) is reported, and an 
approximate relation is given in the 2MASS See- 
ing and Image Statistics ^ to recover the full width 
at half maximum of the point spread function; 

FWHM(arcseconds) = 3.13 sh - 0.46. (2) 

The mean value of this FWHM is about 2.5" for 
the galaxies in our samples (the 2MASS raw im- 
ages have a pixel size of 2" which sets the lower 
limit on the resolution), and the spread is small. 
We derive a correction by numerically convolving 
exponential disks of varying inclination and angu- 
lar sizes with a Gaussian of FWHM = 2.5". Axial 
ratios are measured at a radial distance from the 
galaxy center of 1.7 rd (where refers to the scale 
length of the exponential disk), corresponding to 
r2o in J-band (the ratio r^lriQ is obtained in the 
same way as rd/r2i in Section 5.1). The derived 
correction is parameterized as 

(a/6)corr = (a/6)obs(l - 0.02a; + Q.21x^ - O.Olx^) 

(3) 

with 

X = FWHM ('^\ (4) 

V ^20 / measured 

(where FWHM = 2.5" is used in all cases and 
r2o is measured in arcseconds in the J-band im- 
age). Hereafter, when we refer to the J-band axial 
ratio, we refer to the quantity corrected by the 
above equation. The result of the correction can 
be seen in Figure lb. In a few cases where the 
galaxy is very small and highly inclined the cor- 
rection is large so to avoid over-correcting a max- 
imum value of a/h = 1/0.12 is imposed and these 
galaxies are assumed to have an inclination of 90°. 
From Figure lb it appears that the correction re- 
moves much but perhaps not all of the bias, so 
very small galaxies may appear less inclined than 
they really are. It should also be noted that the 
derived correction assumes light profiles are expo- 
nential, which may not be an ideal approximation 
within the r2o ~ l-7rd radius in J-band where the 
contribution from a bulge may be significant. 

3.3.2. Galactic Extinction and K- corrections 

Magnitudes are corrected for galactic {i.e. orig- 
inating within the Milky Way) extinction using 



fit to many stars in each frame, where /(O) is the 
central surface brightness, r is the radius in arc- 
seconds and a and (3 are free parameters. The 



^2MASS Seeing and Image Statistics (Cutrie 
1998) is an on-line document found at; 
http://spider.ipac.caltech.edu/stciflF/roc/2mass/ 
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the DIRBE galactic extinction maps (Schlcgel, 
Finkbeiner & Davis 1998). The column densi- 
ties listed in the map are converted to extinction 
in the relevant bands by Aj = 0.902E{B - V), 
Ah = 0.576E{B-V) &nd Ak = 0M7E{B-V). A 
small cosmological k-correction is also included us- 
ing a simple linear fit to the low redshift end of the 
models of Poggianti (1997) which is kj = —0.68z, 
kH = -OAOz, and fc^ = -1.522;. 

The magnitude in band X corrected for galactic 
extinction and cosmological k-correction is then 

mx = m^b«-Ax + fcx. (5) 

The isophotal radii also need to be corrected 

for galactic extinction. If we assume the disk pho- 
tometric profiles are exponential with scale length 
rd in band X, then 



(since the surface brightness at radius r in an 
exponential disk is given by /i(r) = /i(0) + 
1.086 r/rd). We also correct sizes for a cosmo- 
logical stretch factor (1 -|- z)'^. 

4. Reddening 

In this section we use the samples defined above 
in Section 2 to investigate the effects of internal 
extinction in the disks by investigating color gra- 
dients with inclination. This test has the advan- 
tage of suffering little from selection effects. For 
the SFI2 sample, we use the I-band axial ratios as 
the measure of inclination since these have errors 
almost three times smaller than the equivalent J- 
band measure from the XSC. The result is shown 
in Figure 2 in which colors (measured within an 
elliptical aperture with a semi-major axis of r2i in 
J-band) are plotted vs. axial ratio. It is clear 
that reddening is present in all colors, which sets 
a lower limit to the extinction in the J and H bands 
by assuming that the extinction in Kg is negligi- 
ble. Following G94 and others, we parameterize 
the internal extinction as 

AM = 7log(a/6). (7) 

The lower limit on 7 is then 0.3 mag in J and 0.1 
mag in H. The formal error in fitting this slope 
to the data is 0.01 mag in both cases, while the 



scatter in the relationship is 0.1 mag. If there is 
extinction present in Kg-band {i.e. 7k ^ 0), the 
actual values of 7 in J and H will be 0.3 -|- 7k 
and 0.1 + 7k respectively. An indication that the 
linear law of Eqn 7 may not be appropriate for the 
highest inclinations is apparent in panels (a) and 
(c) of Figure 2, an issue we'll revisit below and in 
Section 6. 

For the SFI2 sample, I-band magnitudes are 
also available. These magnitudes can be used to 
construct I-J, I-H and I-Kg colors, with the un- 
derstanding that the apertures are not consistent 
between I and the other bands. This mismatch 
should not, in the first order, affect the slope of 
the reddening, just the absolute color. The re- 
sults of this are shown in Figure 3. Under the 
assumption that the extinction in Kg is negligible, 
we derive a lower limit on the extinction in I-band 
of 7 = 0.94 ± 0.03 mag (where the error is the for- 
mal error in determining the slope). G94 find an 
extinction law of 7/ = 1.05 ± 0.08, while Tully et 
al. (1998) derive 7/ ~ 1.0, suggesting that if the 
extinction in Kg is not equal to zero, it is likely to 
be small. 

For the much larger AGC2: sample of galaxies 
(which does not have I-band photometry for all 
the galaxies), we can plot the J-H, J-Kg and H-Kg 
colors against the J-band axial ratios (Figure 4). 
Since the sample is large, we combine the points 
into averages in log(a/6) bins which makes the 
inadequacy of a simple linear law for the relation- 
ship (as noted in Figure 2) stand out. We fit two 
separate lines respectively for log(a/6) less than 
and greater than 0.5, and also include a quadratic 
fit for comparison. The parameterization of AAf 
with a linear law in Eqn 7 is an expedient form, 
and internal extinction can well depart from such a 
linear law, depending on the details of the relative 
distribution of stars and dust (see, e.g. Evans 1992 
and G94). In most sinndations, in fact, AM steep- 
ens with increasing log(a/6) as will be shown for a 
simple photometric model in Section 7. Our adop- 
tion of bi-linear and quadratic fits reflects a least 
departure from simple expediency. The location 
of the break in the bi-linear law at log(a/6) = 0.5 
for all color trends is also practical but arbitrary. 
The slopes of the bi-linear relation are 0.10 ± 
0.02 and 0.31 ± 0.05 for J-H, 0.05 ± 0.02 and 
0.15 ± 0.05 for H-Kg, and 0.15 ± 0.02 and 0.46 
± 0.05 for J-Kg. The quadratic fit is 0.65(20) 
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- 0.03(l)log(a/&) + 0.24(2) [log(a/6)]2 for J-H, 
0.31(25) - 0.04(2)log(a/6) + 0.14(2) [log(a/6)]2 for 
H-Ks, and 0.97(28) - 0.07(2)log(a/6) + 0.39(3) [log(a/6)]2 
for J-Ks. 

5. Effect of inclination on isophotal radius 
5.1. Formal corrections for inclination 

If a galaxy is completely opaque within a given 
isophotal radius in the face-on perspective then no 



No significant dependence of the ratio on inclina- 
tion is observed in any band. We note that the 
f'd/r2i ratios imply that is at only 2.5 ± 0.6 
scale lengths out in the disk, while in H and Kg, 
r2i is still at only 3. 3 ±1.1 scale lengths. The com- 
monly reported isophotal radius at I-band is r23.5 
at 3.5-4 scale lengths which is further out than r2i- 
If the face-on disks were transparent at r2i we 
can rewrite Eqn 9 as 



change in the measured radius is expected with in- 
clination. If it is assumed that the face-on disk is 
transparent at that radius, then the isophotal ra- 
dius should increase as the disk is tilted and we 
look through a longer path length. The expected 
change derived for the simple case of a pure expo- 
nential disk is 



r°,„ + 2.302r21og(a/6) 



(8) 



where r°gg and are respectively the isophotal ra- 
dius and the disk scale length in the face-on case. 
If we assume that the scale length changes with 
inclination like riso (which is justified by the ob- 
servation that the ratio rd/r2i does not change 
with inclination, as seen below) then we derive 

^ = 1 + 2.302— log(a/6). (9) 

We have the scale lengths in I-band measured 
for the SPI2 sample and make the assumption that 
scale lengths do not change significantly between 
I, J, H and Kg-band, as found by de Jong (1996) 
from a sample of 86 face-on spiral galaxies. In 
their sample, scale length is seen to decrease as 
the wavelength increases, but not significantly be- 
tween I and K. A measure of the scale length is 
provided in the XSC for all of the galaxies in our 
sample, but given the notorious difficulties with 
fitting scale lengths, we choose not to use that au- 
tomated measure here, and only use scale lengths 
from the SFI2 sample fit "by hand" (see Haynes 
et al. (1999a) for more details on the SFI2 scale 
lengths). Using these scale lengths we estimate 
rd/r2i = 0.4 ± 0.1, 0.3 ± 0.1 and 0.3 ± 0.1 in J, H 
and Kg for the SFI2 galaxies (where the error here 
is the scatter in the sample), and we assume that 
these ratios hold for all the galaxies in the AGCz 
sample. The ratio rd/r2i in Kg may be slightly 
smaller than in H, but this is within the errors. 



^2iA2i = l+6log{a/b) 



(10) 



where (5 = 0.9 ± 0.3, 0.7 ± 0.2 and 0.7 ± 0.2 in J, H 
and Kg respectively, when we use the rd/r2i values 
given above. The isophotal radius does not change 
{S = 0) if the disk is opaque at r2i so S must 
decrease as the disk becomes less transparent. The 
main limitation on Eqn 9 is that it assumes that 
the photometric profile is exponential, while it is 
expected that the importance of the bulge light 
increases in these NIR bands. 

5.2. Are disks transparent at r2i? 

Here we test the often made assumption that 
the outer disks of galaxies are transparent, espe- 
cially in the NIR. Suppose that stars and dust 
are exponentially distributed with scale lengths 
^d"''* = = rd- The disk opacity at radius r 
and inclination i is then 



T{r,i) = 



r°(0) exp(-r/rd) 



(11) 



where t°(0) is the opacity at r = in the face-on 
perspective. There have been many attempts to 
derive t°(0) for spiral galaxies. G94 derive an up- 
per limit in I-band of Tj (0) < 5 using statistical ar- 
guments with a sample of about 2000 Sc galaxies. 
A lower estimate of t°(0) was obtained by Xilouris 
et al. (1999): they claim that t°(0) < 1 in 7 edge- 
on spiral galaxies they model in K, J, I, V and 
B-bands. Kuchinski et al. (1998) derive values 
of T°(0) using BVRIJHK photometry and radia- 
tive transfer models for 15 spiral galaxies. They 
adopt values in the ranges of 0.15-0.6, 0.1-0.4, and 
0.05-0.2 for J, H and K respectively. Moriondo et 
al. (1998) find t°(0) =0.3-0.5 in H-band from a 
sample of 154 galaxies. Even if Tj(0) were as high 
as 5 it is easy to see that at r23.5,i ~ 3.5r2 the 
disk is close to transparent out to quite high in- 
clinations. However it has been shown above that 
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r2i in J, H and Kg is closer to the center of the 
disk than r23.5,i so (in the above model) even for 
relatively low values of t° (0) the opacity could be- 
come significant in a highly inclined disk. Figure 5 
shows the inclination at which the galaxy becomes 
optically thick (r = 1) plotted as a function of the 
radial distance for r°(0) varying between 0.25 and 
10. 

At r2i in J-band, which is 2.5 ± 0.6 r^, the disk 
can become opaque before it is fully inclined for 
even the smaller r°(0). If we interpret the change 
of slope for high inclinations shown in Figures 4 
and 8a as the onset of r > 1 at r2i, a "break" for 
J-band (r2i,j ~ 2.5rd) and H and Kg (r2i,H/K ^ 
3.3rd) near log(a/6) ~ 0.5 would suggest that the 
disk opacity in the J, H and Kg bands, while low, is 
non-negligible and that even at r2i an important 
fraction of the light is extincted when disks are 
edge-on. 

Recent studies suggest that the dust scale 
length is larger than that of the stars, see for ex- 
ample Moriondo et al. (1998) who use r^^^ = ^r^ 
with ^ = 1.5. In the above model this simply 
changes the position of r2iwith respect to the dust 
scale length so that for example r2i,j ~ 2.5r^ ^ 
(2.5/^)rjJ^''*. In that case a lower central face-on 
opacity is required to explain the behavior. 

5.3. A statistical test for the inclination 
dependence of isophotal radius 

It is possible to test the mean change in isopho- 
tal radius with inclination in a statistical way us- 
ing large samples of galaxies. For the test de- 
scribed below, the biases introduced by the se- 
lection function of the sample can be quite se- 
vere. In order to keep track of these effects, we 
use only galaxies north of declination = —3° from 
the AGCz, since that region is known to be more 
homogeneously and completely sampled in the cat- 
alog. Most of the galaxies in the AGO were orig- 
inally identified by eye from photographic plates 
on which small, dim galaxies are easier to pick out 
if they are edge-on. In an attempt to remove this 
bias we also restrict the sample to B < 15.5. The 
samples are then further reduced to be flux lim- 
ited (using TOoxt from the XSC) or angular size 
limited (using r2i) within each band. The; final 6 
subsamples are taken to represent complete flux or 
angular size limited samples, and are summarized 
in Table 1. 



For either size or flux-limitcd samples, the mean 
value of Tiso = r2i, expressed in kpc, will increase 
with distance. In order to remove this simple man- 
ifestation of the Malmquist bias, wc divide our 
sample into five distance groups with the same 
number of galaxies in each. Within each group the 
mean of r2i is then computed in bins of log(a/6). 
These bins all have an equal number of galaxies 
in them (~30-60 depending on the size of the sub- 
group). Since there arc few galaxies at high in- 
clinations this restricts analysis of subgroups to 
log(a/6) < 0.7-0.8 in most cases. In each of the 
bins the values of < r2i >iog(a/b) a-re scaled by the 
mean of the value for the whole distance group 
(< r2i >dist)- A relation ci -|- C2log(a/6) is fit to 
the points and the value of 5 (sec Eqn 10) is cal- 
culated as C2/C1. The results in Kg-band for both 
the size and flux-limited sub-samples are shown 
in Figure 6, and arc summarized for each band in 
Table 2 along with the values of S derived in Sec- 
tion 5.1 for the case of an exponential disk which 
is transparent at r2i. 

A comparison of panels 6a and 6b, constructed 
from flux and angular size-limited samples respec- 
tively, well illustrates the amplitude of the bias 
introduced by the character of the sample. This 
is explained as follows; as disks become more in- 
clined, on average r2i moves out and the total 
magnitude dims (assuming that the galaxies are 
neither totally transparent throughout the whole 
disk or totally opaque at r2i). In a sample de- 
fined by a threshold angular size, rmin, galaxies 
that in the face-on perspective would be smaller 
than Tmin (and thus excluded from the sample) at 
high inclinations will appear to have a size greater 
than rmin and will be included in the sample. As 
a result, this sample will exhibit a shallow rjgo v. 
log(a/6) slope as shown by Figure 6b. Alterna- 
tively, in a magnitude limited sample, intrinsically 
larger galaxies fall below the threshold magnitude 
if they are edge-on due to the dimming caused by 
extinction. This will have the effect of steepening 
the slope, as seen in Figure 6a. 

Note that if disks were opaque, r2i would not 
change with inclination {6 = 0), and this result 
would be recovered by the analysis of an angu- 
lar size-limited sample. On the other hand, if the 
disk dimming due to internal extinction is small 
or negligible, a flux-limited sample will return the 
correct value of 6 through an analysis such as that 
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Table 1 

Samples for the study of riso and mext vs. log(a/6). 



Band 


Limit 


No. galaxies 


J 


rriext < 11-5 mag 


1774 


J 


r2i > 20 " 


3120 


H 


mcxt < 11 mag 


2274 


H 


r2i > 25 " 


3137 


Ks 


rriext < 10.5 mag 


1733 


Ks 


r2i > 30 " 


2377 



Table 2 

Measured and estimated values of S. 



Band Flux limited Angular size limited Transparent exponential disk 



J 1.2± 0.1 0.5 ± 0.1 0.9 ± 0.3 

H 0.8± 0.1 0.3 ± 0.1 0.7 ± 0.2 

Kb 0.7± 0.1 0.2 ± 0.1 0.7 ± 0.2 



represented in Figure 6a. The presumption of very 
low dimming at Kg-band, supported by the results 
discussed in Section 4, is further supported by the 
similarity of the value of 5 in columns 1 and 3 of 
Table 3 for Kg-band. A final choice on the value 
of S in each band needs to be made after an esti- 
mate of the amount of internal extinction so will 
be discussed in Section 6.3. It is probably safe to 
assume there is very little extinction in Kg so we 
can conclude here that Sk = 0.7 ±0.1. 

6. Effect of inclination on magnitude 
6.1. Formal corrections for inclination 

Note that if the disks of galaxies are completely 
transparent all the way to the center, then there 
should be no change in the total magnitude with 
inclination [i.e. 7 = 0) , while the magnitude 
enclosed within a fixed isophote will brighten as 
the physical radius at which it is measured moves 
out. If the disks arc completely opaque, then both 
the total and isophotal magnitudes will dim by an 
amount 2.5 log(a/6). 



6.2. A statistical test for the inclination 
dependence of total magnitude 

Wo test for the inclination dependence of the 
total magnitude (here taken to be mext as de- 
scribed in Section 3.3) using a method similar to 
that described in Section 5.3. We fit a relation 
AM = 7log(a/6) to the mean values of the mag- 
nitude binned in log{a/b) within five separate dis- 
tance groups. We use the flux and angular size 
limited samples described in Section 5.3, since (as 
for the test in Section 5.3) the interpretation must 
take into account bias introduced by the selection 
function of the sample. For an angular size limited 
sample with threshold rmin, if the disk is not com- 
pletely opaque all the way through, then 7 will be 
increased as intrinsically small and dim galaxies 
move above Tmin at high inclinations. If the disks 
are completely opaque, this type of sample would 
recover the exact result of 7 = 2.5 as the angu- 
lar sizes are unaffected by inclination. For a flux 
limited sample, and if the galaxies are not com- 
pletely transparent, internal extinction will dim 
galaxies at higher inclinations and force 7 to be 
smaller. If the disks are completely transparent, 
this sample would recover the exact result of no 
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change of total magnitude with inchnation. The 
change in 7 with sample selection is observed in 
all photometric bands explored here (as summa- 
rized in Tabic 3), providing further evidence that 
galaxies are neither completely opaque nor com- 
pletely transparent in these bands. We note that 
under the assumption of little extinction in Kg, 
the flux-limited sample in that band should give a 
value of 7 very close to the real value. This implies 
7k = 0.2 ± 0.1. In the other bands, J and H, we 
use this test to provide upper and lower limits to 
the actual extinction law as parameterized by 7. 

6.3. The Modified Holmberg Test 

Following G94 we use a modified version of 
the Holmberg test (Holmberg 1958, 1975). De- 
fine E' = TOiso + 5 log Tiso as the surface brightness 
within an elliptical aperture with semi-major axis 
riso- The original Holmberg test consists of mea- 
suring the change of E' with disk inclination. Un- 
fortunately a similar trend can be found in both 
opaque and transparent galaxies. If the galaxy is 
totally opaque we expect that the radius at which 
a given isophote is measured will not shift with 
inclination and that the magnitude will dim by 
2.5 log(a/&), so a slope of 2.5 is expected. In the 
fully transparent case, fixed brightness isophotes 
will move to larger radii as the galaxies are tilted 
which will also decrease E'. From Eqn. 8 this 
geometric dilution alone {i.e. not taking into ac- 
count the change in the isophotal magnitude as 
the aperture changes) is given by 



51og(r2i,a/r2i,b) = 5 log 



l + 2.302(r>2^)log(«/fe)a' 



l + 2.302(r°/r2°i)log(a/6)b 
(12) 

between axial ratios {a/b)a and {a/b)^. This would 
give a change in E' between log(a/6) = and 1 of 
1.1 ± 0.3 for the rajrix derived in Section 5.1 in 
H and Kg and 1.4 ± 0.3 in J (where the errors in- 
dicate the scatter in rd/r2i). These values provide 
an upper limit to the change since the isophotal 
magnitude will brighten slightly as the aperture 
increases. 

Here, we do not use the isophotal magnitude, 
but rather a measure of the total magnitude (ex- 
trapolated magnitude from the XSC), to remove 
the effects of the physical size of the aperture 
changing with inclination, riso is respectively r2i,j, 
r2i H and r2i k for the three 2MASS bands. In the 



following discussion we will refer to this radius as 
r2i, with the understanding that it is separately 
measured in each of the three bands. The slope fit 
to E' = rrioxt + 51ogr2i vs. log (a/ 6) for galaxies 
in the AGCz sample with lower inclinations (be- 
low log(a/6) = 0.5) is 2.0 ± 0.1 in J, 1.6 ± 0.1 
in H and 1.3 ± 0.1 in Kg. Since this is higher 
than the values given above for the transparent 
case we take it to be indication that there is ex- 
tinction present in the inner parts of these galax- 
ies, while the outer parts are largely transparent. 
This is supported by the radial variation of col- 
ors and surface brightness seen in large galaxies 
in the XSC (Jarrett et al. 2003). If the disk is 
transparent at the edges but not in the center, the 
change of E' = mext + 5 log riso from face-on to 
edge-on will be larger than the values for geomet- 
ric dilution alone and can in fact be larger than 2.5 
as there will be some dimming of the magnitude 
which is dominated by light from the inner parts 
of the disk. 

In order to remove the effect introduced by an 
inclination dependence of the isophotal radius, we 
adopt a modified version of the Holmberg test in 
which the isophotal radius in the face-on view is 
derived via Eqn 10. The value of S which is used in 
Eqn 10 can significantly change the resulting slope. 
This is illustrated in Figure 7 which shows the 
standard Holmberg test (for galaxies observed in 
Kg-band) as well as the test with corrections using 
S = 0.6, 0.7 and 0.8 (spanning the error range for 
5k discussed in Section 5). 

We will discuss the obvious up turn of these re- 
lationships further in Section 7 but for now deal 
with the low inclination (log(a/6) < 0.5) and high 
inclination (log(a/6) > 0.5) regions separately. At 
the low inclinations in Kg-band, 7 changes from 
0.31 (with 6 = 0.6) to 0.05 (6 = 0.8); a sig- 
nificant difference within the plausible range of 
S. As a consequence of this interdependency, the 
best values of S and 7 must be decided upon 
together. Considering all the available evidence 
(Sections 5.1, 5.3 and 6.2) we adopt S = 0.7 in 
Kg-band, 5 — 0.7 in H-band, and S = 0.9 in J- 
band. As a comparison, Moriondo et al. (1998) 
derive Sa = 0.85 ±0.16 at the same isophotal level 
(from their sample of 154 galaxies). The results 
for just these values of 5 in J and H can be seen 
in Figure 8, while the test with the chosen value 
of 6 in Kg-band is shown in Figure 7c. We thus 
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Table 3 
Measured values of 7. 



Band 


Flux limited 


Angular size limited 


J 


0.2± 0.1 


1.2 ± 0.1 


H 


0.1± 0.1 


0.8 ± 0.1 


K. 


0.2± 0.1 


0.7 ± 0.1 



derive 7j = 0.48 ± 0.15, = 0.39 ± 0.15, and 
7k = 0.26 ± 0.15 for log(a/6) < 0.5 from this 

test (where the error indicates the variation in the 
value within the adopted value of i5 ± 0.1. Similar 
values of 7h = 0.16 ±0.22 (Moriondo et al. 1998) 
and 7K 0.22 (Tully et al. 1998) have been found 
before, although Graham (2001) and others sug- 
gest 7i<- ~ 0. While the variation of 7 with the 
wavelength band depends on the detailed form of 
the geometric mixing of stars and dust, it is inter- 
esting to note that, in the simple geometric case 
of a foreground absorbing screen model, the ex- 
tinction curve of our galaxy would yield 7j =0.5, 
7h = 0.3, and 7k = 0.2, if 71 = 1.0 (see G94, Tully 
et al. 1998, Schlcgcl et al. 1998). 

At the higher inclinations (above log(a/6) = 
0.5) the modified Holmberg test is not appropriate 
as the disk is probably not completely transpar- 
ent at r^i- The values of 6 used to correct to the 
face-on value of r2i are therefore too large which 
artificially lowers the change in the average sur- 
face brightness and leads to an underestimate for 
7. The standard Holmberg test {i.e. 5 = 0) does 
not work cither. As discussed above, if there is 
significant extinction present in the inner parts of 
the disks where most of the light comes from, but 
the outer edges are still at least partly transpar- 
ent, the slope of S vs. log(a/6) can be larger than 
2.5 and so is obviously an overestimate of 7. This 
is observed in J, H and Kg (the slopes are 2.9±0.2, 
2.4±0.2 and 2.4±0.2 respectively), providing evi- 
dence that the most inclined galaxies are not to- 
tally opaque or transparent at r2i, and that there 
is significant extinction in the central regions. 
Since the extinction is lowest in Kg we expect that 
the disk is closest to transparent at r2i in this 
band and therefore adopt the slope of the modi- 
fied Holmberg test here to give a tentative value of 
7k — 1.1 ± 0.2 for the more inclined galaxies (up 



to log(a/6) ~ 0.85). We then use the slopes of the 
reddening in J-Kg and H-Kg (above log(a/6) = 0.5: 
7, - 7K = 0.46 ± 0.05 and 7h - 7k = 0.31 ± 0.05 
from Figure 4.) to suggest 7j ~ 1.6±0.2 and 7h — 
1.4 ± 0.2. The difference between the bi-linear 
and quadratic fit is quite small (as can be seen in 
Figure 8), so we also provide a smooth fit to the 
points giving 7j = -0.12(15) -f 1.14(15) log(a/6), 
7h = -0.25(15) 1.12(15) log(a/6) and 7k = 
-0.53(15) -h 1.23(15) log(a/6). Below log (a/6) = 
0.1,0.2, and 0.4 in J, H and Kg respectively the 
value for the correction should be zero - the up- 
turn in the quadratic fit at low inclination is ob- 
viously unphysical. 

7. Photometric Model 

Here we investigate a simple photometric model 
which can be used to explain the upturn in the 
AM vs. log(a/6) relation at high log(a/&) seen in 
all three bands. The model used is that described 
in Section 9 of G94 which is based on the "triple 
exponential" model of Disney et al. (1989). This 
model assumes that the dust is distributed in the 
disk with an exponential scale length equal to the 
stars and that the z distributions are also expo- 
nential with a ratio Q between the scale height of 
the dust and the stars. The face-on central sur- 
face brightness, t°(0) can be specified, as can the 
ratio, q, between the scale height and scale length 
of the stars. We use the "standard galaxy" of G94 
which has 7% of the light in a central bulge. 

Figure 9 shows the difference between the mag- 
nitude at a given axial ratio and the face-on mag- 
nitude calculated for various inputs to the model. 
In the model the upturn becomes more noticeable 
with increasing q {i.e. as the disk becomes thicker) 
and as the central face-on opacity decreases. The 
ratio between the thickness of the disk in the dust 
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and stars has a much smaller effect as shown in 
Figure 9b. 

Figure 10 shows our "best fit" to the Holm- 
berg test points (Section 6.3). Note that at high 
inclinations the points are artificially lowered (see 
discussion in Section 6.3). The simple photometric 
model described here can explain the observed be- 
havior. While the dependence of the models on the 
dust-to-stars thickness ratio is weak, best fit mod- 
els favor low values of (. This is consistent with the 
expectation that stellar disks may appear some- 
what thicker in NIR light, as they are dominated 
by an older stellar population. The models shown 
in Figure 10 were computed for = 0.5. The other 
parameters of the model are t°(0) = 1.1, q = 0.14 
in J-band, r°(0) 0.7, q = 0.14 in H-band and 
r°(0) = 0.3, q = 0.15 in K-band. Reassuringly, 
the central face-on opacity of the best fit model 
decreases with increasing wavelength. Given the 
uncertainties in the Holmberg test (see Section 
6.3) this model should be taken with caution, but 
does provide a plausible explanation for what is 
observed. 

8. Luminosity dependence 

Here we use the AGCz sample of 15,244 galax- 
ies to search for luminosity dependences of the ex- 
tinction laws, as was done in both Giovanelli et 
al. (1995) and Tully et al. (1998). As a mea- 
sure of the total luminosity of the galaxy we use 
Mext.K — 5 log h in K-band since this is least af- 
fected by internal extinction. The luminosity dis- 
tribution is shown in Figure 11. The peak of the 
distribution is at Mext.K — 51og/i = —23, which is 
near to M^^ in K-band (Kochanek et al. 2001) and 
most of the galaxies in the sample are within one 
magnitude of this. 

The expected dependence is that brighter 
galaxies will show more internal extinction at a 
given inclination than dimmer ones. This is due 
both to the increase in physical size (and thus 
optical path lengths) and to the increase in metal- 
licity (and thus presumably dust content) in the 
more luminous galaxies. 

8.1. Isophotal radius 

Here we derive a parametric relation for the 
change of 6 in r2i/r2i = 1 -|- Slog{a/b) with total 
luminosity as measured by Mext.K — 51og/i. We 



discussed the expected value of S for an exponen- 
tial disk with scale length ra which is transparent 
at r2i in Section 5.1. This is given by 2.302rd/r2i 
and is shown by the points in Figure 12 for galaxies 
from the SFI2 sample in 9 bins of luminosity (each 
with an equal number of galaxies). There is a clear 
trend for the ratio rd/r2i (and thus 6 if the disks 
are transparent at r2i) to decrease with luminos- 
ity, as was also found by Giovanelli et al. (1995) 
in I-band, who also discuss why this trend is to 
be expected. We fit a cubic to these points which 
we truncate at both the highest and lowest lumi- 
nosities. These fits are shown by the solid lines 
in Figure 12 and the fit parameters are listed in 
Table 4. 

8.2. Magnitude 

Here we apply the "modified Holmberg test" 
discussed in Section 6.3 to the full AGC2; sample 
to investigate the dependence of internal extinc- 
tion (as parameterized by 7 in Eqn 7) on galaxy 
luminosity. We divide the sample into 10 luminos- 
ity bins with ^ 1500 galaxies per bin. We correct 
the isophotal radius to the face-one view using the 
cubic fit to the points in Figure 12. We then fit 
a single linear law to E vs. \og{o./h) for each lu- 
minosity bin as the subdivision in luminosity bins 
restricts statistics, especially at high inclinations. 

The extent of the luminosity dependence of 7 
depends critically on the value S used to correct 
radii to the face-on values. In the modified Holm- 
berg test, an over-estimate of S leads to an under- 
estimate of 7 (and vice- versa). Since 5 clearly de- 
creases with luminosity, using a constant value of 6 
for all galaxies thus leads to an over-estimate of 7 
at low luminosities and an under-estimate at high 
luminosities and suppresses the expected luminos- 
ity dependence. The cubic fit in Figure 12 follows 
luminosity dependence of 6 (under the assumption 
of transparency at r2i) and should be used here. 
The values of 7 derived by using this S are shown 
in Figure 13 as a function of luminosity. We find 
that over much of the range of luminosities the 
value of 7 is roughly constant in each band. In 
support of this we observe that the slope of the 
reddening does not change much with luminosity. 
It is clear that at low luminosities the extinction 
in these bands does decrease significantly and (un- 
like as found for galaxies in I-band by Giovanelli et 
al. 1995) does fall off to zero for galaxies less lumi- 
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nous than Mext,K - 51og/i = -20, -20.7 and -20.9 
in J, H and Kg respectively. The sohd Hne in Fig- 
ure 13 shows our adopted luminosity dependence 
for 7, which is broken into two regions. The rela- 
tion is 77 = -0.26a; + 0.77, 7// = -0.26a; + 0.58, 
and 7k = -0.13x + 0.31, for x > 0.8, and jj = 
-0.04X+0.60, 7ff = -0.01x-H0.38, and -/k = 0.20, 
for X < 0.8 (where x = Mext,K -5\ogh + 23). It 
should be noted that a straight line can fit within 
the errors on the points (dotted line in Figure 13), 
which are ±0.15 in 7. 

If the assumption of transparency at r2i is 
not valid then we are over-estimating 6 and thus 
under-estimating 7. High luminosity galaxies have 
more extinction, but r2i/rd is also larger in these 
galaxies which could mean that they are more 
likely to be transparent at r2i than low luminos- 
ity galaxies. What bias the assumption of trans- 
parency might introduce into this determinations 
of the luminosity dependence of 7 is not clear. 

9. Conclusions 

Two samples of spiral galaxies are constructed, 

both of which have photometry from 2MASS and 
redshifts from the AGC. For the smaller sample 
(which has over 3000 galaxies), I-band photome- 
try is also available. We apply a variety of sta- 
tistical tests to look for disk inchnation depen- 
dences in these galaxies. The results of these tests 
can be highly sensitive to the selection function of 
the sample and some interdependencies exist, but 
taken together they provide reliable information 
on the average internal extinction in the galaxies. 
We parameterize internal extinction via Eqn 7 and 
change of the isophotal radius via Eqn 10; these 
results are summarized in Tables 4 and 5, where 
the numbers between brackets indicate the error 
on the last digit of the parameter value. 

From the trends of NIR color with inclination in 
the SFI2 sample, we derive a firm lower limit to the 
total extinction in J and H-band of 7j = 0.3 + 7k 
mag and 7h = 0.1 + 7k mag respectively. The 
scatter in these trends is on the order of 0.1 mag, 
while the statistical error on determining the slope 
is 0.01 mag. Using the same argument, we derive a 
lower limit on the extinction in I-band of 0.9 -|- 7k 
mag, which when compared with those derived by 
G94 and TuUy et al. (1998), is consistent with 
very low extinction in Kg-band. Using the much 



larger AGCz sample, wc note that the simple lin- 
ear relation given by Eqn 7 does not adequately 
describe the variation of the 2MASS colors with in- 
clination, but an upturn at high inclinations takes 
place instead which we fit with a bi-linear law. 

The radius r2i is seen to be closer to the center 
of the disk in J, H and Kg (~ 2-3 scale lengths) 
than r23.5 in I-band (at 3.5-4 scale lengths), and 
we work through a simple model to estimate at 
what inclination the disk might become opaque 
at that radius for a given central face-on opacity. 
Using samples defined according to different se- 
lection criteria, we explore the inclination depen- 
dence of r2i and compare the data with the depen- 
dence that would be expected for pure exponential 
disks which are transparent at r2i at all inclina- 
tions. We find that complete transparency at r2i is 
inconsistent with the data of very nearly edge-on 
galaxies. Upper and lower limits are estimated for 
6 (as defined by Eqn 10), and we adopt as most 
probable values 5j ~ 0.9 ± 0.1, Sn = 0.7 ± 0.1 and 
Sk = 0.7 ±0.1. 

Tracking the flux dependence on inclination and 
using a modified version of the Holmberg test we 
investigate the amount of internal extinction in 
disks in the three bands. This analysis confirms 
that a simple value of 7 (as defined by Eqn 7) is 
not adequate to describe extinction for all values 
of log(a/fe). Because of this wc fit a bi-linear law 
fitting separately high and low inclination systems 
as separated by log(a/6) = 0.5. The value of 7 
calculated from this test is sensitive to how the 
isophotal radii are corrected to the face-on value 
{i.e. the value of S that is used in Eqn 10) and 
can change by a significant amount. Using the 
adopted values of 6 we find 7j = 0.48 ±0.15, 7h = 
0.39±0.15, and 7k = 0.26±0.15 for log(a/6) < 0.5. 
If 6 is increased, 7 decreases and vice- versa. 

At the higher inclinations the modified Holm- 
berg test is not appropriate as we expect a much 
smaller variation of r2i with inclination than is 
corrected for. This implies that the measured 
slope is only a lower limit to the value of 7 in 
Eqn 7. We expect that the limit will be closest 
to the actual value of 7 in Kg where the extinc- 
tion is lowest and therefore adopt 7k = 1.1±0.2 + 
0.13/ log(a/6). Using information on reddening we 
then derive 7j = 1.6 ± 0.2 ± 0.24/ log(a/6) and 
7h = 1.4 ± 0.2 + 0.20/ log{a/b). 

We explore a simple photometric model for ex- 
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Table 4 

Summary of results for 5 in r2i/r2i = 1 + ^log(a/6). 



Method J-band H-band Ks-band 

For transparent exponential disks 0.9(3) 0.7(2) 0.7(2) 

Statistical test 

lower limit 0.6(1) 0.4(1) 0.2(1) 

upper limit 1.2(1) 0.8(1) 0.7(1) 

Adopted 0.9(1) 0.7(1) 0.7(1) 

Luminosity dependence 0.87 + O.llx 0.70 + 0.07a; 0.67 + 0.10a; 

(a; = Mext,K-51og/i + 23) +0.05a;^ + 0.007a;^ +0.05a;^ + 0.009a;^ +0.05a;^ - 0.004a;^ 



Table 5 

Summary of results for 7 in AM = 7log(a/6) 



Method 


J-band 


H-band 


Kg-band 


Reddening test 








All SFI2 


0.3(1)+7K 


0.1(1)+7K 






0.2(1)+7H 






AGCz 








low inclination'' 


0.15(2)+7K 


0.05(2)+7K 






0.10(2)+7h 






high inclination* 


0.46(5)+ 0.08/ log(a/6) 


0.15(5) + 0.025/ log(a/6) 






0.31(5) + 0.05/log(a/6) 






Statistical test 








lower limit 


0.2(1) 


0.1(1) 


0.2(1) 


upper limit 


1.2(1) 


0.8(1) 


0.7(1) 


Holmborg tost'' 








low inclination* 


0.48(15) 


0.39(15) 


0.26(15) 


high inclination* 


> 1.6 + 0.24/ log(a/6) 


> 1.3 + 0.20/ log(o/6) 


> 1.1 + 0.13/ log(o/6) 


Adopted 








low inclination'' 


0.48(15) 


0.39(15) 


0.26(15) 


high inclination* 


[1.6(2) + 0.24/ log(a/6)] 


[1.4(2)+0.20/log(a/6)] 


[1.1(2) + 0.13/log(a/6)] 


Luminosity dependence 








[X = Mext.K - 5 log + 23) 








X > 0.8 


-0.26a; + 0.77 


-0.26a; + 0.58 


-0.13a; + 0.31 


X < 0.8 


-0.04a; + 0.60 


-O.Ola; + 0.38 


0.20 



*Limit between low and high inclinations is log(o/6) = 0.50. 
''Assumes 6 =0.9, 0.7 and 0.7 in J, H and Kg respectively. 
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tinction due to dust in spiral galaxies. In the 
model the dust is distributed exponentially both 
radially and in the ^;-direction. The scale length 
of the dust is assumed to be the same as the 
stars, while a ratio ( is assumed between the scale 
heights. The "best fit" models to the data (all 
with C = 0.5) suggest a central face-on opacity, 
T°(0) = 1.1 and intrinsic axial ratio, q = 0.14 
for J-band, t°(0) = 0.7 and g = 0.14 for H-band 
and T°(0) = 0.3 and q = 0.15 for Ks-band. This 
model explains well the upturn in the extinction 
at high inclinations, which becomes more obvious 
for lower values of t°(0). 

In Tables 4 and 5 we give estimates of our fa- 
vored solutions for 5 in Eqn 10 and 7 in Eqn 
7 as compromises of derivations from different 
methods. They are listed in the row labelled 
"Adopted". 

Finally, we explore the luminosity dependence 

of 5 and 7. Under the assumption of transparency 
at r2i the value of 5 decreases with luminosity. 
We use the modified Holmberg test to explore the 
luminosity dependence of 7 and find that the ex- 
tinction drops off to zero below a luminosity of 
Mext.K - 51og/i = -20, -20.7 and -20.9 in J, H 
and Kg respectively, and is roughly constant at 
7 = 0.61, 0.39 and 0.19 at the higher luminosities. 
There is marginal evidence of further increase at 
the very highest luminosities. The adopted lumi- 
nosity dependences are also listed in Tables 4 and 
5. 

We have explored the commonly held view that 
there is little internal extinction in external galax- 
ies at NIR wavelengths, especially Kg-band. We 
have shown that except for the least luminous 
galaxies, there is evidence for internal extinction 
(which increases with luminosity), and that in the 
most inclined galaxies the extinction becomes sig- 
nificant. 

We are very grateful to Dr. Tom Jarrctt for 
his help in accessing the 2MASS Full release Ex- 
tended Source Catalog and for his reading and 
comments on the paper. This publication makes 
use of data products from the Two Micron All 
Sky Survey, which is a joint project of the Univer- 
sity of Massachusetts and the Infrared Processing 
and Analysis Center/California Institute of Tech- 
nology, funded by the National Aeronautics and 
Space Administration and the National Science 



Foundation. This work has been partially sup- 
ported by the NSF Grant No. AST-0098526. 
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Fig. 1. — Isophotal angular semi-major axis (in arcsec) at the 20th magarcsec"^ level plotted against the 
log of the axial ratio (which is measured at about the 20th mag arcsec"^ level), (a) As observed; the line 
shows the theoretical limit from our simulations, (b) Corrected for seeing with FWHM= 2.5" as described 
in text. The obvious banding is due to the accuracy of tabulated values in the XSC. 
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Fig. 2. — NIR color plotted against the log of the axial ratio (measured from I-band images) for 3035 galaxies 
for which I-band photometry is available (SFI2 sample). The slopes of the linear fits are 0.2 in J-H (a), 0.1 
in H-Kg (b), and 0.3 in J-Ks (c). The formal error in determining the slope is 0.01 mag in all cases, while 
the scatter in the relationship is roughly 0.1 mag. Assuming there is no extinction in Kg-band (which will 
be tested below) this gives a lower limit to the value of 7 (where AM = 7log(a/&)) in J and H bands of 0.3 
and 0.1 mag respectively. Note the indication of non-linearity in this relationship (especially in J-Kg (c)) 
which is much more evident in Figure 4 and is discussed in Sections 4, 6.3 and 7. 
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Fig. 3. — NIR color plotted against the log of the axial ratio (measured from the I-band images) for 3035 
galaxies for which I-band photometry is available (SFI2 sample). Note that the aperture in which the I-band 
magnitude is measured is not the same as the J, H and Kg apertures, but this should not affect the slope 
only the offset. The lines shown are fits to the median in bins of log (a/6) which is also shown in plot (d). 
The slope, 7 is 0.73 in I- J (a), 0.82 in I-H (b) and 0.94 in I-Kg (c). The scatter is 0.4 mag in all cases, and 
the error on determining the slope from the medians is 0.03 mag. Assuming no extinction in Kg-band this 
then gives a lower limit to the extinction in I of 0.9 magnitudes, consistent with the value of 7 = 1.05 ± 0.08 
derived in G94 and 7 ~ 1.0 in TuUy et al. (1998). For comparison with Fig 2 the vertical scale is exactly 
three times larger. 
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J-band log(a/b) I— band log(a/b) 

Fig. 4. — NIR color plotted against the log of the axial ratio (measured from the J-band images) for 15224 
galaxies (AGC^r sample) binned in groups of 300 (error bars show the statistical error on determining the 
mean). The fit to the points is given as a quadratic (solid line) and as a bi- linear relation with a break at 
log(a/6) = 0.5 (dotted line). Panel (d) shows log(a/6) in J vs. log(a/6) in I (both corrected for seeing) for 
the 3035 SF2 galaxies in the sample; the scatter is 0.1 in log (a/ 6) and there is no significant indication of 
bias in the J-band axial ratios. For comparison with Figure 2 the range in ordinate value in (a), (b) and (c) 
is exactly three times smaller. 
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Fig. 5. — The observed axial ratio (assuming an intrinsic axial ratio, q = 0.13 for the dotted lines or 0.2 for 
the solid lines, where cos^ i = {{h/af' — — q^)) at which a galaxy with a pure exponential disk of dust 

(with r^"*** = r*i) becomes opaque is shown as a function of galactic radius for various central optical depths 
(0.25, 0.5, 1, 2, 3, 4, 5, 10 from left to right). Far from the center, the disk is transparent at most inclinations, 
but further in the disk can become optically thick for fairly modest r° (0) at high enough inclinations. 
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Fig. 6. — Averaged values of the isophotal radius, r2i in Kg-band in bins of log(a/6) scaled by the value of all 
the galaxies in the same distance bin (as described in Section 5.3). Fit to this is a function ci + C2 log(a/5), 
the plot shows the data and the fit divided by Ci. The slope, S = C2/C1 is 0.7 ± 0.1 in the magnitude limited 
sample (1733 galaxies) and 0.2 ± 0.1 in the radius limited sample (2377 galaxies), where the error is the 
formal error on determining the slope. 
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Fig. 7. — Holmberg test in Kg-band (see Section 6.3). The values of 5 shown are that used to correct the 
isophotal radius to the face-one value using r2i/r2i = 1 + S\og{a/b). A bi-lincar fit is shown (dotted line) 
with the break at log(a/6) = 0.5. The fitted slopes are (a) 1.34 and 2.38, (b) 0.31 and 1.34, (c) 0.26 and 
1.06, and (d) 0.05 and 1.07. The error on the slope is on the order of 0.1 for the low inclinations and 0.2 for 
the high inclinations. Also shown is a quadratic fit for comparison (solid line). This smooth function may be 
more desirable, but the unphysical upturn at low log(a/6) may make the bi-linear law more useful. Points 
are averages in bins of \og{a/b) with roughly 300 galaxies per bin, where the error bars show the error in 
determining the mean in the bin. 
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Fig. 8. — Modified Holmberg test in J and H-bands with the chosen value of S (see Section 6.3). Shown is a 
bi-linear fit to both cases where tlie break is at log(a/b) = 0.5 (dotted line). The slopes arc 0.48 ± 0.1 and 
1.57 ± 0.2 in J-band (a), and 0.39 ± 0.1 and 1.30 ± 0.2 in H-band (b). Also shown (long dashed line) is the 
adopted correction. The solid lines show a quadratic fit for comparison. 
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Fig. 9. — Difference between the "observed" magnitude and the face-on value from the "triple exponential" 
model of Disney et al. (1989) which assumes dust and stars are exponentially distributed with the same 
scale length. We vary the intrinsic axial ratio q, the ratio ( between the scale heights of the dust and stars 
and the central face-on opacity, r°(0). (a) Results for ( = 0.5, q = 0.15 and r°(0) = 0.1,0.5,1.0,1.5,2,4,10 
(from right to left), (b) Results for t°(0) = 0.5 and q = 0.1,0.15,0.2 (right to left). The small effect of 
varying ( between 0.1 (dotted line), 0.4 (short dashed), 0.7 (long dashed) and 1.0 (solid line) is also shown. 
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Fig. 10. — Shown here are the data points from the Holmberg test in J-band (from Figure 8a) in (a), H-band 
(from Fig 8b) in (b) and Kg-band (from Fig 7c) in (c), offset to zero at log (a/6) = 0. The solid lines show 
the photometric model which best fits (by eye) the data. All models are for C = 0.5. For the J-band data 
(a), T°(0) = 1.1 and q ^ 0.14; for H-band (b), t°(0) = 0.7 and q = 0.14 and for K^-band (c), t°(0) = 0.3 
and q = 0.15. 
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Fig. 11. — Distribution of total magnitude in Kg-band of galaxies in the AGC^; sample. 
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Fig. 12. — The change of 2.302rd/r2i with luminosity for galaxies in the SF12 sample. This parameter 
approximates S if disk are nearly transparent at 7'2i. Data points show the averages in 9 bins of luminosity 
(with equal numbers of galaxies), while the solid lines show a cubic fit to the points truncated for the highest 
and lowest luminosity galaxies to be the value of the fit at Mext.K — 51og/i = —21.4 and -24.2 respectively. 
The fits are respectively for Moxt.K — 5 log h > —21.4, between -21.4 and -24.2 and brighter than -24: S = 1.19, 
d = 0.87 + 0.11a; + O.OSa;^ + O.OOlx^, 6 = 0.80 for J-band (a); S = 0.97, 6 = 0.70 + 0.07a: + 0.05.t2 + 0.009a;3, 
5 = 0.66 for H-band (b); and S = 0.94, 5 = 0.67 + O.lOa; + O.OSa;^ - 0.004a;^ 5 = 0.63 for K,-band (c). 
X = Mext.K — 5 log /i + 23 in all cases. 
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Fig. 13. — Values of 7 as a function of luminosity in (a) J-band, (b) H-band and (c) K-band derived using 
the modified Holmbcrg test with 5 as given by the truncated cubic fit to the points in Figure 12. The dotted 
line shows a straight line fit to the points, the solid line our adopted luminosity dependence. 
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